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Abstract 



We analyze a fiber-optic component which could find multiple uses in 
novel information-processing systems utilizing squeezed states of light. Our 
approach is based on the phenomenon of photon-number squeezing of soli- 
ton noise after the soliton has propagated through a nonlinear optical fiber. 
Applications of this component in optical networks for quantum computation 
and quantum cryptography are discussed. 
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Information-processing systems where information is carried by nonclassical states of 
light (e.g., photon number states or squeezed states), though not as yet implemented 
with standard telecommunications in mind, nevertheless present an attractive alternative 
for such novel applications as quantum computation || or secure quantum-key distribution 
||. As envisioned currently such systems can be constructed using linear optics, linear 
mixing of creation and annihilation operators (linear Bogoliubov transformations which 
include squeezing as a special case), and nonlinear operations for state preparation and 
detection. 

One major obstacle to current practical implementation of such systems is the lack of 
components which would preserve the noise characteristics of the transmitted information. 
This requirement is crucial because squeezed states are highly sensitive to loss, as a simple 
analysis reveals. In optical communication networks, optical amplifiers are inserted along 
the fiber links in order to restore exponentially decaying signal power to acceptable levels. 
These devices, represented by an idealized model of a phase-insensitive amplifier (PIA), 
add at least 3 dB of noise with coherent-state inputs when the amplifier gain is high || 
(this minimal noise figure is referred as the Standard Quantum Limit, or SQL [[]]) and thus 
may destroy precisely engineered noise statistics of the transmitted states. The quantum- 
key distribution scheme of Gottesman and Preskill fU, for example, calls for at least 2.51 
dB of squeezing in channels with weak noise, so introducing a PIA into the underlying 
communications infrastructure may severely affect the performance of the scheme. 

In this Letter we concentrate on a component of such information-pocessing systems 
that could prove useful for generating and manipulating squeezed states of light, photon- 
number squeezed (PNS) states in particular. As recent experimental reports indicate 0, it 
is possible to generate sub-Poissonian light (a near-PNS state) by means of the following 
simple setup. A soliton pulse is launched into an optical fiber and then frequency-filtered 
at the fiber exit. The resulting quantum- mechanical state of light exhibits photon- number 
fluctuations below the coherent-state level, with maximum squeezing observed when the 
fiber is three soliton periods long. To illustrate the utility of this setup to our goal, we 
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analyze the arrangement in which the above process is iterated by launching the filtered 
soliton through another fiber-and-filter stage. From here on, we will refer to our setup as 
the Dual-Stage Squeezer (DSS), as opposed to the Single-Stage Squeezer (SSS) described at 
the beginning of the paragraph. 

Our motivation to suggest the use of the DSS in optical networks for quantum informa- 
tion processing comes from the observation, further elaborated below, that the DSS can be 
thought of as the original SSS operating on squeezed (rather than simply coherent) input 
states. This property of the DSS makes it a natural choice for easily implementable sources 
of squeezed states, whose degree of squeezing can be controlled by tuning the soliton-pulse 
parameters, and also for devices that enhance the degree of squeezing introduced into the 
quantum network. For reasons mentioned above, such squeezing enhancers would conceiv- 
ably play a key role in practical implementations of the Gottesman-Preskill quantum key 
distribution scheme using squeezed states. 

In a recent study, Levandovsky et al. || have employed the soliton perturbation approach 
H to obtain a complete theoretical description of the quantum-noise statistics of spectrally 
filtered solitons. This linearization approach is valid whenever the photon-number noise is 
small compared to the average number of photons in the soliton, which is usually the case 
in most experiments. We brifely recount the main idea of their approach and then apply it 
to our analysis of the DSS. 

Nonlinear evolution of an electromagnetic pulse propagating through a lossless optical 
fiber is governed by the quantum nonlinear Schrodinger equation 
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where a(r, £) is the annihilation operator of the field and (£, r) are the dimensionless space 
and time coordinates. The corresponding classical equation has a fundamental soliton solu- 
tion Oo(t, £) = e l ^ 2 sechr = fo(r)e 1 ^ 2 , given here in canonical form with two photons per 
pulse. We write the annihilation operator as 

a(r,0 = [/o(r)+Aa(r,0]e <e/2 ) (2) 



where Aa is the annihilation operator that represents the perturbation of the soliton 
mean field by quantum noise and satisfies the usual equal-space commutation relations, 
[Aa(r,0,Aa(r',0] = [Aat(r, 0, AaV, 0] = 0, [Aa(r, 0, AaV, 0] = S(r - r 1 ), every- 
where inside the fiber. We make the linearizaiton approximation by substituting Eq. (0) 
into Eq. ([I]) and discarding all terms that are 0(Aa 2 ), thus separating the problem into the 

classical NLSE for the mean field ao(r,£) and the linearized operator equation 
f)Ah i r) 2 

= 2 h Ab + 2 ^o(r, \ 2 Ab + ial(r, (3) 
for Ab = Aae 1 ^ 2 . In what follows, we shall disregard phase factors of the form e* 6 ^, — oo < 



9 < oo, because the filtered light is directly detected in a SSS. 

The solution of Eq. (§) can be written as an eigenfunction expansion || 

Aa(r,0 = ^- f[V c (n,0m,r) + V s (n,0fs(n,T)]dn+ £ V5(0/i(r) (4) 

i=n,p,r,& 

with operator coefficients, where the discrete eigenmodes /„, f p , f T , and fg represent per- 
turbations of the soliton mean field due to changes in photon number, momentum, time, 
and phase respectively; and f c and f s are the symmetric and anti-symmetric continuum 
eigenmodes that represent perturbation of the dispersive radiation in the fiber. Detailed 
analysis of these modes, along with their time-domain and frequency-domain forms, can be 
found in Ref. 0. 

The ^-dependent Hermitian operators Vt, i G M = {c, s,n,p,r,9}, are obtained by 
projecting Eq. (|) onto the eigenmodes {/, | i G M} of the equation adjoint to Eq. @ which 
is obtained by reversing the sign of the Ab^ term. The relevant orthogonality relations are 
(fii fj) = Ajj, where the inner product is defined by (/, g) = Re / f{r)g*{T)dr\ Ay = 5ij in 
all cases except for i = j G {c, s}, where A^- = 2ir5(fl — Q'). 

Defining the time-domain cosine quadrature operator Aa c = (Ao + Ao')/2, we write the 
time-domain correlation function G(r, r';£) = 4(Aa c (r, ()Aa c (r', £)) and the corresponding 
covariance function C(r, r';^) = G(r, t';£) — 4(Aa c (r, ())(Aa c (r', ^)). Assuming that the 
filter i^(o;) at the fiber exit is linear and imposing the realizability condition < < 1, 

we obtain the perturbation operator after the filter from the frequency-domain relation 



Aa ont (u,0 = \H(u)\Aa(u,0 + y/l- \H(u)\H(u), (5) 

where v is a vacuum-state operator associated with the frequency-dependent loss due to the 
filter. The observed squeezing S(£) is quantified by normalizing the output photon-number 
variance to the average output photon number: 
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where Cn(uj, lo'\ £) = C(u,w';£) — 2tt5(uj — uj') is the normally ordered part of the Fourier 
transform C(u;,u/;£) of the covariance function C(r, t';£), and /o( w ) = 7rsech^ is the 
Fourier transform of fo{r). 

As stated by Levandovsky e£ a/. ||, in the case of the SSS the correlation function 
G(u,lu';£) is equal to the covariance function C(uj, because the perturbation at the 
fiber entrance (£ = 0) is white coherent-state quantum noise. The observed squeezing in 
dB, given by —10 log is shown in Fig. [IJ for the case of a bandlimited parabolic filter 
|l0l H(u) = 1 — uj 2 /r] 2 , \u\ < rj, where the bandwidth r\ is adjusted to give 10% loss. The 



analytical expression for the correlation function G in time domain, as well as the observed 
squeezing for other types of filters, are given in Ref. [§]. 

Now we turn to the analysis of the DSS. In order to render the problem tractable while re- 
taining its essential physical features, we have made the following assumptions: (a) the filter 
in the first stage is weak enough, so that the quantum-mechanical average of its output may 
be treated as the soliton mean field ao(r, £) plus a small perturbation A6o(r, £) = Aa (r)e ! ^ 2 ; 



(b) all relevant frequency- domain quantities are narrowband [IT], so that the number of pho- 
tons at high frequencies that are cut off with a bandlimited filter is negligible, and the mean 
field can still be treated as a fundamental soliton; and (c) as the pulse propagates through 
the second stage, the quantum-mechanical average (Aa(r)) remains independent of £ and 
is equal to Aa (r). This last assumption can be given precise mathematical meaning by 
ensuring that the maximum steady-state (r — ► oo) error which results from assuming that 
Aa (r)e 1 ^ 2 is a solution of the linearized Eq. ([|) is vanishingly small fL2l . With these as- 
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sumptions in place, the DSS becomes equivalent to the SSS, but with a squeezed-state input 
and an additional filter inserted between the soliton source and the entrance to the fiber. 

Since the filter is assumed to be weak, we may write its transfer function in the form 
H(u) = 1 + h(u), where h(u) is small, so that Aa (w) = h(u)fo(uj) is the frequency-domain 
form of the quantum-mechanical average (Ad(r)) which, as we have assumed, is independent 
of £. Therefore, at the fiber entrance to the second stage, the soliton mean field is perturbed 
by an ideal squeezed state with the quantum- mechanical average Aao(r) and the squeezing 
parameter r determined by the fiber length and the filter transfer function of the first stage. 
It is then easy to see that the time-domain covariance function for the DSS is given by 

Cdss(t,t';0 = e- 2r G S ss(r,r';0 - 4Aa (r)Aa (r'), (7) 

where Gsss( r ; T '] 6) is the correlation function for the SSS and £ = at the fiber entrance to 
the second stage. The observed squeezing Sbss(0 a ^ the output of the DSS is then computed 
by substituting the normally ordered covariance function Cdss,at(^, w'; £) into Eq. (|6]). 

We have analyzed numerically the following arrangement. If the fiber in the first stage is 
three soliton periods long, then with our particular filter the maximum squeezing is 2.8 dB, 
which corresponds to r ~ 0.32. The initial perturbation Aciq(u)) is given by — ^foyoj) (for 
all frequencies \uj\ < f], a, condition which is sufficient in accordance with our narrowband 
assumption). The maximally squeezed output of the SSS is then launched through the 
second stage. The observed squeezing at the output of the DSS with the first stage generating 
maximally squeezed output is shown in Fig. as a function of the fiber length in the second 
stage, with maximum squeezing of 6.1 dB observed if the fiber in the second stage is three 
solition periods long. 

The results of the preceding analysis suggest that by iterating the SSS it is possible to 
construct a wide variety of components for generation and manipulation of squeezed states 
in optical networks for quantum communication and computation. As already mentioned, 
the SSS can be used to generate near-PNS states whose statistical parameters can be tuned 
by varying the pulse power and width. The main incentive to use the DSS for generation of 
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near-PNS states is the simplicity of implementation: in order to get more squeezing, we just 
add another stage to the SSS with another parabolic filter, thus avoding the use of optimized 
filters whose shapes are complicated and difficult to implement. |8j In other words, compared 
to a SSS, one gets more squeezing out of a DSS while using less fiber and easily realizable 
filters. 

As for manipulation of squeezed states in such networks, the DSS, with more than a 
twofold increase in squeezing (as measured in decibels), renders an excellent illustration of 
how one can use fiber nonlinearity to enhance the noise statistics of light transmitted through 
the network. Such an enhancement would be needed, e.g. if the squeezing introduced pre- 
viously was degraded by loss. To summarize, the capabilities offered by fiber nonlinearity 
need to be closely explored with such novel applications as quantum information processing 
in mind. 

The authors acknowledge useful discussions with D. Levandovsky. This work was supported 
in part by the U.S. Army Research Office through MURI grant DAAD19-00-1-0177. 
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where /(r, £) is the approximate solution of the c-number version of Eq. (§) and £[■} is 
the linear operator acting on the (mean) perturbation on the right-hand side of Eq. @. 
In our case, /(r,£) = (I/77 2 ) exp(i£/2)<i 2 sech r/<iT 2 , and the error is given by 

6 sech 3 r tanh 2 r 

^steady = hm 5 = 0. 

In fact, the quantitiy under the limit sign rolls off very rapidly in r. 



^P-£[/(r,0] 
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FIGURES 

FIG. 1. Observed squeezing (in dB) vs. fiber length (in soliton periods) in the SSS. The filter 
has a bandlimited parabolic frequency response giving 10% loss. 

FIG. 2. Observed squeezing (in dB) vs. second-stage fiber length (in soliton periods) in the 
DSS with the first stage generating maximally squeezed output (the first-stage fiber is 3 soliton 
periods long). The filter has a bandlimited parabolic frequency response giving 10% loss. 
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Fig. [3] - "Generation and manipulation ..." by Raginsky and Kumar. 
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